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We present an update of a previous global analysis of the experimental data on azimuthal asymmetries in semi- 
inclusive deep inelastic scattering (SIDIS), from the HERMES and COMPASS Collaborations, and in e + e~ — ► 
h\h2X processes, from the Belle Collaboration. Compared to the first extraction, a more precise determination 
of the Collins fragmentation function and the transversity distribution function for u and d quarks is obtained. 



1. Introduction 

The study of the nucleon spin structure has 
recently made remarkable progress. Our under- 
standing of the longitudinal degrees of freedom, 
concerning both the intrinsic motion and the spin 
content of partons inside unpolarized and longi- 
tudinally polarized fast moving nucleons, respec- 
tively encoded in the unpolarized quark distribu- 
tion, g(cc), and the helicity distribution, Aq(x), is 
now quite accurate. 

A full knowledge of the nucleon quark struc- 
ture in the collinear, k± integrated, configura- 
tion cannot be reached without information on 
the third twist-two parton distribution function: 
the transversity distribution, ATq(x) (also de- 
noted as hi{x)). Despite its fundamental impor- 
tance [T] and the intense theoretical work of the 
last decade [2j, this function has only very re- 
cently been accessed experimentally. The main 
difficulty in measuring transversity is that, being 
a chiral-odd quantity, it decouples from inclusive 
deep inelastic scattering (DIS), since perturbative 
QED and QCD interactions cannot flip the chi- 
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rality of quarks. 

The only way to access this distribution is by 
coupling it to another chiral-odd quantity. To 
such a purpose one can look for a chiral-odd part- 
ner either in the initial or the final state. In the 
first case the most promising approach, and the 
cleanest one from the theoretical point of view, 
is the study of the double transverse spin asym- 
metry, Att, in Drell-Yan processes. This mea- 
surement is in principle feasible at RHIC, but the 
small x region covered at a cm. energy y/s — 200 
GeV and the fact that in pp collisions one mea- 
sures the product of two transversity distribu- 
tions, one for a quark and one for an antiquark, 
lead to Att values of the order of a few per- 
cents [3]. A much larger A T t, around 20-40%, 
could be observed in Drell-Yan processes in pp 
interactions at s ~ 200 GeV 2 , as proposed by 
the PAX Collaboration [4151617] . However, this 
requires the availability of polarized antiprotons, 
which is an interesting, but formidable task in it- 
self. Other double transverse spin asymmetries 
for inclusive production of photons or pions are 
strongly suppressed by the large gluon contribu- 
tion in the unpolarized cross sections [8|9j . 

For the case of a chiral-odd partner in the final 
state, we mention the spin transfer in p^p — > X 
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or £p^ — > ^'A^AT processes, where the final hy- 
peron acts as a polarimeter. Here A^g couples 
to another unknown quantity: the fragmentation 
function for a transversely polarized quark into a 
transversely polarized baryon. 

Recently a promising suggestion has been 
made and is going to be pursued: the combined 
study of two-hadron inclusive production in single 
polarized DIS and in electron-positron annihila- 
tion measurements |10) . In this context the new 
quantity which appears is the so-called dihadron 
fragmentation function describing the hadroniza- 
tion of a quark in two hadrons [ll|12|13|14|15j . 

Meanwhile, the most accessible and fruitful 
channel is the azimuthal asymmetry A^^ h+< ^ s ^ 
in SIDIS processes, namely lp^ — > IhX, involv- 
ing the convolution of the transversity distribu- 
tion with the Collins fragmentation function [16] : 
the spin and transverse momentum dependent 
(TMD) function parameterizing a left-right asym- 
metry in the fragmentation of a transversely po- 
larized quark into an unpolarized hadron. 

This study has been and still is under active 
investigation by the HERMES, COMPASS and 
JLab Collaborations. 

A crucial breakthrough has been achieved 
thanks to the independent measurement of the 
Collins function (or rather, of the convolution of 
two Collins functions), in e + e~ — > h\hiX un- 
polarized processes by the Belle Collaboration at 
KEK [17]. By combining the SIDIS data from 
HERMES QJ] and COMPASS [15], with the Belle 
data, a global fit leading to the first extraction of 
the transversity distribution and the Collins frag- 
mentation functions was performed in Ref. [20] . 

Recently, much higher statistics data on these 
spin azimuthal asymmetries have become avail- 
able: the HERMES Collaboration have presented 
charged and neutral pion, as well as kaon az- 
imuthal asymmetries [3T]; the COMPASS Col- 
laboration have presented their measurements, 
still on a deuteron target, but now for sep- 
arate charged pion and kaon production |22j : 
the Belle Collaboration have issued new high- 
precision data of the Collins asymmetry in e + e~ 
annihilation [23j . 

Therefore, we reconsider here our previous 



analysis and study the impact of these new data 
on the extraction of both the transversity and the 
Collins functions. 

2. Formalism 

We recall here the main steps necessary to cal- 
culate the azimuthal asymmetries in SIDIS and 
in e + e~ annihilation, addressing to Ref. [20] for 
details. This approach is based on an extension of 
the ordinary collinear factorization theorems with 
inclusion of a new class of spin and TMD distribu- 
tions 24 25 26 27J. For such processes k± factor- 
ization has been proven |28|29|30] in the regime 
of low observed transverse momenta (compared 
to the large scale of the processes). Another im- 
portant result, crucial for this analysis, is the uni- 
versality of the Collins function entering SIDIS 
and e + e~ processes, as discussed in Refs. [31|32| . 
We will restrict ourselves to tree level expressions, 
as currently done in phcnomenological studies, 
neglecting the soft factor coming from gluon re- 
summation [28)29130] responsible of potential Su- 
dakov suppression [33] . We will briefly comment 
on this in the sequel. 

2.1. SIDIS 

Let us consider the single spin asymmetry for 
the process ip 1 ' — > i'hX: 

(ftptp 1 ->e'hx _ ( j6 a e P l ^e'hx 
da^ — da^ 

do-T + da-i ' ^' 

where d 6 a ip ' ~ * thx = da^^ is a shorthand nota- 
tion for {d 6 a epUl ^ ehx )/(dxdydzd 2 P T d(l)s) with 
x,y,z the usual SIDIS variables. The 4>s depen- 
dence originates from the cross section depen- 
dence on the angle between the proton (trans- 
verse) polarization vector and the leptonic plane. 
We adopt here the standard SIDIS kinematics 
according to the "Trento Conventions" [34], see 
also [20] ■ By considering the sm(4>h + 4>s) mo- 
ment of Ajjt, w e are able to single out the ef- 
fect coming from the spin dependent part of the 
fragmentation function of a transversely polar- 
ized quark (embedded in the Collins function, 
A N D h / q r(z,p±) or H^~ q (z,p±) [33]) coupled to 
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the TMD transversity distribution (Axq(x,k±))- 
More explicitly, we get 



A 



'>-{<t>h+<t>s) 



UT 



= 2 



J d<j>s d(f>h [d<J T - da 1 ] sin(^, l + 4> s ) 
J d(f> s d<t>h [dal + da±] 
E 9 eg ^rq(x, fcj_) ® A N D h/qT (z,p ± ) 

Eg e l fq/p( X , k ±) ® D h/q (z,P±) 



, (2) 



where <£> stands for a convolution on the trans- 
verse momenta (see Eq. (4) of Ref. [20] for full 
details and related comments). 

This analysis can be further simplified by work- 
ing at 0(k±/Q) and adopting a Gaussian and 
factorized parameterization of TMDs. In particu- 
lar for the unpolarized parton distribution (PDF) 
and fragmentation (FF) functions we use: 



f q /p(x,kx) = f q /p(x) 



Dh/ q (z,p±) = D h/q (z) 



-kl/(kl) 

,-p1/(p1) 



(3) 
(4) 



with (kj_) and fixed to the values found in 
Ref. [35] by analyzing unpolarized SIDIS: 

(k 2 x ) = 0.25 GeV 2 , (pi) = 0.20 GeV 2 . (5) 

Integrated parton distribution and fragmenta- 
tion functions, f q / p (x) and D h / q (z), are avail- 
able in the literature; in particular, we use the 
GRV98LO PDF set [36] and the DSS fragmenta- 
tion function set [37] , 

For the transversity distribution, Axq(x,k±), 
and the Collins FF, A N D h / q t(z,p±), we adopt 
the following parameterizations [2"U] : 

Arqix,^) = i/V* (x)[f q/p (x) + Aq(x)} 



(6) 



7r(fc2) T 



A N D h/q r(z,p ± ) = 2Af q c (z)D h/q (z) 

e -pl/(pl 



xh(px) 



(7) 



with 



a ( a + B)( a +0) 
M? (x) = N? x a (l - xf 



h + s) 



(7+<5) 



(9) 



h{p 1 _) = V¥e^-e-Al M K (10) 

and -1 < Nf < 1, -1 < N$ < 1. We assume 
(kj) T = (k\). The helicity distributions Aq(x) 
are taken from Ref. [38]. Notice that with these 
choices both the transversity and the Collins func- 
tion automatically obey their proper positivity 
bounds. 

Using these parameterizations we obtain the 
following expression for A s ^^ h+< ^ s ' ': 

a ut =—T—m — G (^) 



M h 1 + (1 - yf 



E 9 e 2 A^(x) + A^z)] Af q c (z)D h/q (z) 



Y. q e lfq/p(. X ) D h/ q {z) 



(11) 



where C(Pt) is given by 



^- (» 2 i ) 2 p- P t/( P t) c (P 2 ) 



with 



M 2 h ( P 2 



,(PtU ={pl\ C] +z 2 (kl) 



(13) 



When data or phenomenological information at 
different Q 2 values are considered, we take into 
account, at leading order (LO), the QCD evo- 
lution of the transversity distribution. For the 
h/qt, as its scale dependence is 



Collins FF, A N D 



unknown, we tentatively assume the same Q 2 evo- 
lution as for the unpolarized FF, D h / q . 

By performing a best fit of the measurements of 
HERMES, COMPASS and Belle Collaborations 
we then fix the free parameters, a, (3, 7, 8, N q , N q 
and Mh appearing in A^y^'* + ^ s ' ) (q = u,d). 

2.2. e + e~ — > h-ih^X processes 

One might think that hadron production in 
e + e~ collisions is the cleanest process for the 
study of TMD polarized fragmentation functions, 
like the Collins function, thanks to the lack of cor- 
responding TMD effects in the initial state. How- 
ever, in the process e + e~ — > qq there is no trans- 
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verse polarization transfer to a single, on-shell fi- 
nal quark. Therefore, the single Collins effect, 
i.e. the asymmetry in the distribution around 
the jet thrust axis (given by the fragmenting 
quark direction) of hadrons produced in the quark 
fragmentation, cannot be measured. Instead, in 



hadron production from e + e 



2 jets 



events, the Collins effect can be observed when 
the quark and the antiquark are considered simul- 
taneously. The Belle Collaboration at the KEK- 
B asymmetric-energy e + e~ storage rings have 
in fact performed a measurement of azimuthal 
hadron-hadron correlations for inclusive charged 
dihadron production, e + e~ — ► ir + ir~X [17|23j . 
This asymmetry has been interpreted as a direct 
measure of the Collins effect, involving the con- 
volution of two Collins functions. 

Two methods have been adopted in the exper- 
imental analysis performed by Belle. These can 
be schematically described as (for details and def- 
initions see, e.g., Refs. |39I20I23| ): 
i) the "cos(^i method" in the Collins-Soper 
frame where the jet thrust axis is used as the z 
direction and the e + e~ — > qq scattering defines 
the xz plane; 

it) the "cos(2y>o) method", using the Gottfried- 
Jackson cm. frame where one of the produced 
hadrons (J12) identifies the z direction and the 
xz plane is determined by the lepton and the h 2 
directions. There will then be another relevant 
plane, determined by z and the direction of the 
other observed hadron hi, at an angle ipo with 
respect to the xz plane. 

In both cases one integrates over the magni- 
tude of the intrinsic transverse momenta of the 
hadrons with respect to the fragmenting quarks. 
For the cos(^i + ^2) method the cross section for 
the process e+e~ — > h\h 2 X reads: 

^.e+e - ->h 1 h 2 X 



dz\ dz2 d cos 9 d(ip± 
3a 2 
4s 



P2) 



^7 ^2 e l {(! + cos 2 9) D hl/q ( Zl ) D h2 /q(z 2 ) 



+ - 



sin 2 9 



cos^i+^a) A N D hl/q r( Zl ) A N D h2/ql (z 2 )} 1 

(14) 

where 9 is the angle between the lepton direction 



and the thrust axis and 

A N D h/q r(z) = J d 2 P± A N D h/q i(z, P± ) 



(15) 



Normalizing to the azimuthal averaged unpo- 
larized cross section one has: 



A 12 {z 1 ,z 2 ,9, <pi + (p 2 ) 

\ ^p-e+e - -^hihzX 



(da) dz\ dz 2 dcos9d(ipi + ip 2 ) 
„2 a 



1 

= 1 + T 



4 1 + cos 2 9 



cos(^i + ip 2 ) 



E g e 2 A^ l/9T (zQA^ 2 /^2 



(16) 



For the cos(2<y9o) method, where the Gaussian 
ansatz becomes extremely helpful, the ana- 
logue of Eq. (fl6|) reads 



A (z 1 ,z 2 ,9 2 ,tpo) 
= 1 + 1 2122 



sin 2 2 



n z( + zj 1 + cos 2 9 2 



cos(2<p ) 



Z q e 2 q A N D hl/q r(z 1 )A N D h3/9 r(z 2 ) 

T l q e2 q D h 1 /q(zi) D h2/q (z 2 ) 



(17) 



where 9 2 is now the angle between the lepton and 
the h 2 hadron directions. 

To eliminate false asymmetries, the Belle Col- 
laboration [23] consider the ratio of unlike-sign to 
like-sign pion pair production, Ajj and A^. 

For fitting purposes, it is usually convenient to 
express these relations in terms of favoured and 
unfavoured fragmentation functions, 



Ar+/-M 



D 
D 



fav • 



/ S,S 



(18) 
(19) 



and similarly for the A D's. 

3. Results 

A combined fit of SIDIS asymmetries together 
with e+e~ -> h x h 2 X data, Eqs. (|11|16I17|) . allows 
the simultaneous extraction of the transversity 
distribution and the Collins fragmentation func- 
tions. We assume flavour independent values of 
a and (3 (neglecting transversity distributions of 
sea quarks) and, similarly, we assume that 7 and 8 
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are the same for favoured and unfavoured Collins 
fragmentation functions; we then remain with a 
total number of 9 parameters. 

The first study along this line was presented in 
Ref. 20J . Here we repeat the analysis, exploit- 
ing the new high-precision data recently released 
by the HERMES [2TJ and COMPASS [22] Col- 
laborations for SIDIS, and by the Belle Collab- 
oration [23] for e + e~ annihilation processes, in 
order to refine and reduce the uncertainty of the 
previous extraction. 
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Figure 1. Fit of HERMES [21] data. The shaded 
area corresponds to the statistical uncertainty in 
the parameter values, see text. 



New data from COMPASS operating on a 
transversely polarized hydrogen target have re- 
cently been released [ID]: these are not included 
in the fit but compared with our predictions. 

The two sets of Belle data, coming from two 
analyses of the same experimental events, are not 
independent. Therefore we include only one set 
of data in the fit, either A$ or A\2 data. In this 
analysis we report the results obtained by using 
A\2 data, the cos(<^i + ip 2 ) method. The conse- 
quences of fitting Aq instead of A12 are presently 
under investigation. 

In Figs. [1] and [2] we show the best fit to the 
HERMES [2] and COMPASS [22] data, respec- 



Figure 2. Fit of COMPASS [22] data. The shaded 
area corresponds to the statistical uncertainty in 
the parameter values, see text. The extra tt phase 
in addition to 4>h + 0s comes from the different 
convention adopted by COMPASS. 



tively. Notice that the tt° data (HERMES) have 
not been used in the fit; in Fig. [T] we show our 
estimates, based on the extracted transversity 
and Collins functions, and compare them to data. 
Fig. [3] shows the fit to the Belle A\ 2 asymmetry, 
whereas in Fig. [5] our predictions for the A asym- 
metry are compared with data [23 . 

The curves shown are evaluated using the cen- 
tral values of the parameters in Table [T] while 
the shaded areas correspond to a two-sigma devi- 
ation at 95.45% Confidence Level (for details see 
Appendix A of Ref. [41]). 

Tabled] collects the results of our best fit to the 
new data sets [21122123] . while in Figs. O and [SJ we 
show our updated transversity distribution and 
Collins fragmentation functions together with the 
uncertainty bands of our previous extraction [20] . 
We can definitely say that the two extractions are 
compatible with each other, with the new error 
bands strongly reduced. The transversity for up 
quarks results now larger (compared to our previ- 
ous extraction) , while that for down quarks is bet- 
ter constrained in sign and non compatible with 
zero. In this respect the new data from SIDIS 
have been crucial. It is worth noticing that while 
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Table 1 

Best values of the free parameters for the u and 
d transversity distribution functions and for the 
favoured and unfavoured Collins fragmentation 
functions. We obtain x 2 /d.o.f. = 1.3. No- 
tice that the errors generated by MINUIT are 
strongly correlated, and should not be taken at 
face value. The significant fluctuations in our re- 
sults are shown by the shaded areas in the plots. 
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Figure 3. Fit of the Belle [23] data on the A 12 
asymmetry (the cos(y>i + ^2) method). 



the transversity for up quarks is strongly con- 
strained by HERMES data, in particular through 
the positive pion azimuthal asymmetry, the addi- 
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Figure 4. Comparison of our predictions with 
Belle [23] data for the A§ Belle asymmetry (the 
cos(2<po) method). 



tion of COMPASS deuteron data to the fit allows 
a better determination of Axd. We recall here 
that, in analyzing SIDIS data, we have assumed 
the transversity distributions for sea quarks and 
antiquarks to vanish. The extracted Collins FFs 
are well constrained and much smaller than their 
positivity bounds, with the unfavoured Collins 
function large in size and negative, consistently 
with other extractions [42I43I20] . 

A word of caution has to be added here since 
SIDIS data (HERMES and COMPASS) are col- 
lected at a much smaller scale (Q 2 ~ 2.5 GeV 2 ) 
compared to the Belle data (Q 2 = 110 GeV 2 ). 

Both azimuthal asymmetries in SIDIS and in 
e + e~ collisions involve spin and TMD functions 
whose behaviour upon scale variation should be 
described in the context of Collins-Soper factor- 
ization [28.30J. Beyond tree level this would result 
in a soft factor entering TMD convolutions, with 
the corresponding Sudakov suppression. This, as 
discussed in Refs. |44|45j . might imply an under- 
estimation of the Collins function as extracted at 
tree level from the azimuthal asymmetry at Belle. 
Hence the combined extraction of the transversity 
from SIDIS at a lower Q 2 (less Sudakov suppres- 
sion), might lead to an overestimation of A^g. 
This issue is currently under study. Here, as in 
Ref. 20], the Q 2 dependence of the Collins FF 
is included assuming it to be the same as that 
of the unpolarized fragmentation function, D h / q : 



Update on transversity and Collins functions from SIDIS and e + e data 



7 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



x kJGeV) 

Figure 5. The transversity distribution functions 
for u and d flavours as determined by our global 
fit, at Q 2 = 2.4 GeV 2 ; we also show the Sof- 
fer bound [46] (highest or lowest lines) and the 
(wider) uncertainty bands of our previous extrac- 
tion [2U] . 
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z p.(GeV) 



Figure 6. Favoured and unfavoured Collins frag- 
mentation functions as determined by our global 
fit, at Q 2 — 2.4 GeV 2 ; we also show the positiv- 
ity bound and the (wider) uncertainty bands as 
obtained in Ref. [20] , 



although this might not be the proper evolution, 
it should mitigate the above-mentioned effect. 

As it is well known, in a non rclativistic the- 
ory the hclicity and the transversity distributions 
should be equal. We then show in Fig. [7] the 
extracted transversity distribution together with 
the helicity distribution of Ref. [3H] at Q 2 = 2.4 
GeV 2 . It results that, both for u and d quarks, 
we have \Axq\ < \Aq\. 

Another interesting quantity, related to the 
first x-moment of the transversity distribution, 
is the tensor charge: 

5q= f dx (A T q - A T q) = f dx A T q (20) 
Jo Jo 

where the last equality is valid for zero antiquark 
transversity, as assumed in our approach. From 
our analysis we get, at Q 2 = 0.8 GeV 2 , 

Su = 0.54lg;g Sd = -0.23±g;^ . (21) 

Such values are quite close to various model pre- 
dictions 47 48 49 50 for tensor charges which 
span the ranges 0.5 < Su < 1.5 and —0.5 < 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



x kJGeV) 

Figure 7. Comparison of the extracted transver- 
sity (solid line) with the helicity distribution 
(dashed line) at Q 2 = 2.4 GeV 2 . The Soffer 
bound [46] (blue solid line) is also shown. 
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our result 
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Figure 8. Tensor charge from different models 
compared to our result. 1: Quark-diquark model 
of Ref. [47], 2: Chiral quark soliton model of 
Rcf. 05J, 3: Lattice QCD 05], 4: QCD sum 
rules [50] , 

5<i < 0.5 (see Fig. |5J. In this context it is worth 
mentioning a subtle point concerning the strong 
scale dependence of the tensor charge, recently 
addressed in Ref. [51] . For the effective models 
of baryons, as those referred to above, the choice 
of their starting energy scale and their Q 2 evolu- 
tion could play a significant role and, eventually, 
mask the true nature of the model. Consequently, 
the results shown in Fig. [8] where our LO phe- 
nomcnological extraction seems in better agree- 
ment with the quark-diquark model of Ref. |47] 
than with other models, should be taken with 
some care. A safer quantity, totally scale indepen- 
dent, and therefore easy to compare with, would 
be the ratio of two tensor charges. From our fit, 
for instance, we obtain 8d/8u = — 0.42+q'2o° 3 , and 
all model predictions considered above would fall 
within our uncertainty band, as shown in Fig. 7 
of Ref. [5]]. 

4. Predictions 

We now use the extracted transversity and 
Collins functions to give predictions for new mea- 




0.5 1 1.5 

P T (GeV) 



Figure 9. Predictions for the single spin asymme- 
try v 4™( ( ^+^s+7r) com p arec [ t Q preliminary data 
by the COMPASS experiment operating with a 
transversely polarized hydrogen target 00] , 



surements performed or planned at COMPASS 
and JLab. The transverse single spin asymme- 
try A^^ h+ ^ s ^ has been recently measured by 
the COMPASS experiment operating with a po- 
larized hydrogen target (rather than a deuterium 
one). In Fig.[9jwe show our predictions compared 
with these preliminary data. The agreement is 
excellent. 

In Fig. [TU] we present our estimates for JLab 
operating with a proton target at 12 GeV. Notice 
that JLab results will give important information 
on the large x region, which is left basically un- 
constrained by the present SIDIS data from HER- 
MES and COMPASS. In this region our estimates 
must be taken with some care. We recall that 
the large x behaviour of our parameterization is 
controlled by the same f3 parameter for Atu and 
Ard (since present data do not cover the large 
x region). The same is true for the Collins frag- 
mentation functions, whose large z behaviour is 
driven by the same parameter S for favoured and 
unfavoured Collins FFs. On the other hand for 
the small to medium x region, well constrained by 
SIDIS measurements, data support the choice of 
a universal behaviour x a for Atu and Axd. The 
future JLab measurements, which will extend to 
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Figure 10. Estimates of the single spin asymme- 
try A^:^ h+ ^ s ^ for JLab operating with proton 
target. 



The extracted transversity distributions and 
the Collins fragmentation functions allow to com- 
pute the azimuthal asymmetry A S j^^ h+ ^ s ^ for 
any SIDIS process. In particular, our predictions 
for the COMPASS measurements with a proton 
target are in very good agreement with prelimi- 
nary data, while the large x behaviour of A^q, 
yet unknown, could be explored by JLab experi- 
ments. These will provide further important tests 
of our complete understanding of the partonic 
properties which are at the origin of SSAs. 

Further expected data from Belle will allow to 
study in detail not only the z dependence of the 
Collins functions, but also their p± dependence. 

The TMD approach to azimuthal asymmetries 
in SIDIS and e + e~ — ► h\hi X processes has defi- 
nitely opened a new powerful way of studying the 
nucleon structure and fundamental QCD proper- 
ties. 



larger x values, will test the validity of this ap- 
proximations. 

5. Conclusions 

We have performed a re-analysis of recent high- 
precision experimental data on spin azimuthal 
asymmetries which involve the transversity dis- 
tributions of u and d quarks and the Collins frag- 
mentation functions. The values of the 9 free pa- 
rameters are fixed by simultaneously best fitting 
the HERMES, COMPASS and Belle data. 

All data can be accurately described, lead- 
ing to the extraction of the favoured and un- 
favoured Collins functions, in agreement with 
similar results previously obtained in the litera- 
ture [43 42 20 . In addition, we have improved the 
extraction of the so far poorly known transver- 
sity distributions for u and d quarks, Atu(x) and 
Axd(x). They turn out to be opposite in sign, 
with |Atc?(2:)| smaller than |Arw(cc)|, and both 
smaller than their Soffer bound [46] . The pre- 
vious uncertainty bands are strongly reduced by 
the present analysis. The new distributions are 
compatible with our previous extraction |20] and 
close to some model predictions for the transver- 
sity distribution. 



Acknowledgments 

U.D. would like to thank the organizers of this 
very interesting workshop for their kind hospital- 
ity. We acknowledge the support of the European 
Community - Research Infrastructure Activity 
under the FP6 "Structuring the European Re- 
search Area" program (HadronPhysics, RII3-CT- 
2004-506078). M.A., M.B., and A.P. acknowledge 
partial support by MIUR under PRIN 2006 and 
by the Helmholtz Association through funds pro- 
vided to the virtual institute "Spin and strong 
QCD" (VH-VI-231). 

REFERENCES 

1. J. P. Ralston and D.E. Soper, Nucl. Phys. 
B152 (1979) 109. 

2. V. Barone, A. Drago and P.G. Ratcliffc, 
Phys. Rep. 359 (2002) 1, |hep-ph/0104283| 

3. O. Mar tin et al., Phy s. Rev. D60 (1999) 
117502, |hep~ph/9902250| 

4. P AX, V. Ba rone et al., (2005), 
|hep-ex/0505054 

5. M. Anselmino et a l., Phys. Lett. B594 (2004) 
97, hep-ph/0403lT4| 

6. A.V. Efremov, K. Goeke and P. 



10 



M. Anselmino et al. 



Schw eitzer, Eur. Phys. J. C35 (2004) 
207, jhep-ph/0403124] 

7. B. Pasquini, M. Pincetti and S. Boffi, Ph ys. 
Rev. D76 (2007) 034020, [hep^ph/0612094 

8. A. Mukherjee, M. Stratmann and W. Vo- 
gelsang, Phys Rev. D67 (2003) 114006, 
|hep-ph/0303226| 

9. A. Mukherjee, M. Stratmann and W. Vo- 
gelsang, Phys Rev. D72 (2005) 034011, 
|hep-ph/0506315| 

10. A. Bacchetta et al., (2008). [arXiv:0812.061ll 
[hep-ph]. 

11. A.V. Efremov, L. Mankiewicz and N.A. Torn- 
qvist, Phys. Lett. B284 (1992) 394. 

12. J .C. Collins an d G.A. Ladinsky, (1994), 
|hep-ph/9411444| 

X. Artru and J.C. Collins Z. Phys. C69 
(1996) 277, |hep-ph/9504220| 
R.L. Jaffe, X.m. Jin and J. Tang, Phys. Rev. 
Lett. 80 (1998) 1166, |hep-ph79 709322. 
A. Bacchetta and M. Radici, Phys. Rev. D69 
(2004) 074026, |hep^ph/0311173[ 
J.C. Collins, Nucl. Phys. B396 (1993) 161, 



13. 



14. 



15. 



16. 



17. 



18. 



19. 



20. 



21 



22. 



23. 



24. 



25. 



26. 



27. 



28. 



hep-ph/9208213. 

Belle, K. Abe et al., Phys. R ev. Lett. 96 
(2006) 232002, [hep^ex/0507063[ 
HERMES, A. Airapetian et al., Phys. 



Rev. 



Lett. 94 (2005) 012002, |hep-ex/ 0408013. 

COMPASS, E.S . Ageev et al., Nucl. Phys. 

B765 (2007) 31, |hep~ex/0610068| 

M. Ans elmino et al., Ph ys. Rev. D75 (2007) 

054032, |hep-ph/0701006| 

HERMES, M. Diefenthaler, (2007), 

larXiv:0706.2242l [hep-ex] . 

COMPASS, M. Alekseev et al., (2008), 

larXiv:0802.2i"60l [hep-ex] . 

Belle, R. Seidl et al., Phys. Rev. D78 (2008) 

032011. larXiv:0805.2975l [hep-ex], 

P. Mulders and R. Tangerman, Nucl. Phys. 

B461 (1996) 197, |hep-ph/9510301| 

D. Boer and P.J. Mulders, P hys. Rev. D57 

(1998) 5780, |hep-ph/9711485| 

A. Bacchetta et al, JHEP 02 (2007) 093, 

|hep-ph/0611265| 

U. D'Alesio and F. Murgia, Prog. Part. Nucl. 
Phys. 61 (2008) 394, arXiv:0712.4328l [hep- 
ph]. 

J.C. Collins and D.E. Soper, Nucl. Phys. 



B193 (1981) 381. 

29. X.-d. Ji, J.-p. Ma and F. Yuan, P hys. Lett. 
B597 (2004) 299, |hep-ph/0405085| 

30. X.-d. Ji, J.-p. Ma and F. Yuan, Ph ys. Rev. 
D71 (2005) 034005, [nep^ph/0404183[ 

31. A. Metz, Phys. Lett. B549 (2002) 139. 

32. J.C. Collins and A. Metz, Phys. R ev. Lett. 
93 (2004) 252001, [hep^ph/04082491 

33. D. Boer, Nuc l. Phys. B603 (2001) 195, 
|hep-ph/0102071| 

34. A. Bacchetta et al., Phys. 
117504, |hep-ph/0410050| 

35. M. Anselmino et al., Phys 
074006, |hep-ph/0501196| 

36. M. Gluck, E. Re ya and A. Vogt 
J. C5 (1998) 461,| hep-ph/9806404| 

37. D. de Florian, R. Sassot and M. Strat- 
mann, Phys Rev. D75 (2007) 114010, 
|hep-ph/0703242| 

38. M. Gluck et al., Phys. Rev. D63 (2001) 



Rev. D70 (2004) 
Rev. D71 (2005) 
Eur. Phys. 



094005, hep-ph/0011215 

39. D. Boer, R. Jakob and P.J. Mulders, N ucl. 
Phys. B504 (1997) 345, |hep-ph/9702281| 

40. COMPASS, S. Levorato, (2008), 

larXiv:0808.0086l [hep-ex]. 

41. M. Anselmino et al, (2008). larXiv:0805.2677l 
[hep-ph] . 

42. A.V. Efremov, K. Goeke and P. 
Schweitzer. Phys. Rev. D73 (2006) 094025, 



hep-ph/0603054 



43. W. Vogelsang and F. Yuan, Phys. Rev. D72 
(2005) 054028, [hep^ph/0507266[ 

44. D. Boer, Nu cl. Phys. B806 (2009) 23, 
larXiv:0804.2408l [hep-ph]. 

45. D. Boer, (200S). larXrv:0808. 2886 [hep-ph]. 

46. J . Soffer, Phys Rev. Lett. 74 (1995) 1292, 
|hep-ph/9409254| 

47. I.C. Cloet, W. Bentz and A.W. Thomas, 
Phys. Lett. B659 (2008) 214. larXiv:0708.3246l 
[hep-ph]. 

48. M. Wakamatsu, Phys. Lett. B653 (2007) 398, 
larXiv:0705.29T7l [hep-ph]. 

49. QCDSF, M Gockeler et al. , Phys. Lett. B627 
(2005) 113, |hep-lat/0507001| 

50. H.-x. He and X.-d. Ji, Phys. Rev. D52 (1995) 



2960, |hep-ph/94 12235 

51. M. Wakamatsu, (2008), larXiv:08 11.41961 
[hep-ph] . 



